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CSJ . Using the NRGR effective field theory formalism we calculate the remaining source multi- 

2jJ' pole moments necessary to obtain the spin contributions to the gravitational wave amplitude 

.^ ' to 2.5 Post-Newtonian (PN) order. We also reproduce the tail contribution to the waveform 

linear in spin at 2.5PN arising from the nonlinear interaction between the current quadrupole 



and the mass monopole. 



INTRODUCTION 



Direct detection of gravitational waves (GWs) is an international experimental effort which 
includes the earth-based LIGO /Virgo detectors l|] and the future space-based mission eLISA [4]. 



(N 
> 

a^ 

. ' Among the most promising sources for these detectors are coalescing compact binary systems, 

m ' 

<^ • whose dynamics can be studied using Post-Newtonian methods for the inspiral phase and numerical 

(N 

simulations for the merger. 

The search for GWs from compact binaries is performed by the method of matched filtering, 
^ I which relies on having accurate theoretical templates for the expected GW signal. Matched filtering 

for the periodic GW signals emitted by binary systems is most sensitive to the phase. That is the 
reason for the common use of restricted waveforms, which include only the leading PN effects in the 
amplitude but higher PN orders in the phase. However, a precise knowledge of the GW amplitude 
has several benefits: R increases the mass range through the inclusions of higher harmonics 



J.U.C 



improves parameter estimation [JJ, yields a higher angular resolution and a better the measurement 
of the luminosity distance ^ , and is necessary for the comparison of numerical and analytic results 
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Measurements of black hole spins indicate that black holes in binary systems may frequently 



be close to maximally rotating 



3]. Thus it is timely and necessary to develop more accurate GW 



templates that include spin effects to higher orders in both the phase and the amplitude. For 



nonrotating compact objects, the GW amplitude from binary inspirals is known to 3PN 
spin effects have only been fully included at 1.5PN order and partially at 2PN [9[]. 



d], but 



In this paper we compute the spin effects in the GW amplitude to 2.5PN. We use the NRGR 



effective field theory (EFT) formalism developed in jlCll4l6l| (see [ITI. Il8l| for reviews), with which 
the conservative part of the dynamics for rapidly rotating compact objects has been completed 



to 3PN order in 
have been a 
reported in 



19l427l|. (The next-to- leading order (NLO) spin-orbit and spin-spin Hamiltonians 



so computed in [28l43Cll| using more traditional PN methods, and full agreement was 
31[. The conservative spin(l)-spin(2) potential was recently computed at NNLO using 



both the EFT 



21 and the ADM formahsm (33|.) 



Blanchet, Buonanno and Faye computed in [34] the source multipole moments linear in spin 
for the spin effects in the energy flux and the GW phase to 2.5PN in a different framework and 
with different conventions from ours. In 35| we calculated the source multipole moments needed to 
obtain the spin contributions to the power emitted in GWs to 3PN order. These results include the 
corrections linear in spin, previously obtained in [3j], and the hitherto uncalculated spin(l)-spin(l) 
and spin(l)-spin(2) contributions. The latter provided the last ingredients required to determine 
the GW phase evolution to 3PN order, which will be presented in [3a] where we will include a 
detailed comparison to the results of ^34] . 



In the present paper we extend our work in 



351 ] and calculate the missing ingredients to obtain 



35| 



the spin effects in the gravitational waveform to 2.5PN order, which were not computed in 
since they enter the GW phase beyond 3PN. The source multipole moments calculated here are 
linear in spin, and consist of the LO current and mass 16-poles, LO current 32-pole and the NLO 
current octupole. The former three multipole moments were previously computed in 3j] which 
derived the LO multipoles linear in spin to all orders in the multipole expansion, whereas the NLO 
current octupole linear in spin is derived here for the first time. See Table jTI] below for a clear 
display of which m ultip ole moments are needed for the waveform to 2.5PN. We also reproduce the 



recently computed 37| spin contribution from the tail effect at 2.5PN. 



To 



22 



gether with numerical simulations 



24 



, the results reported here and elsewhere 



11 



13 



IS- 



351 ] provide the building blocks for the construction of more accurate GW templates to be 



used in the data analysis of rapidly rotating compact binaries, e.g. 



mUm- 



This paper is organized as follows. In the next section we briefly review the main ingredients for 
our computation of the spin components in the GW amplitude within the EFT formalism. In the 
ensuing section we provide the results for the remaining spin induced source multipole moments, 
which contribute to the waveform to 2.5PN order. Finally, we discuss nonlinear gravitational 
effects in the GW amplitude and reproduce the tail effect due to the interaction between the 
current quadrupole and the mass monopole, which is linear in spin and also enters at 2.5PN 
order. A detailed discussion of the formalism to calculate nonlinear effects in the waveform will be 
elaborated upon in 42l |. 



We will rely heavily on previous work and thus encourage the reader to consult the literature 
on the subject for more details IOtIJ, Il9l422l. |35|. Throughout the paper we work in ^ = c = 1 
units, with Newton's constant and the Planck mass related by G = 1 / {?>2T[mi) . 

2. GRAVITATIONAL WAVEFORM, MULTIPOLE MOMENTS AND SPIN IN NRGR 

2.1. Gravitational waveform 

The gravitational waveform has a complicated structure due to the nonlinear nature of general 
relativity. However the leading effects in the PN expansion are given by the linear propagation of 
GWs emitted by the source multipole moments. The waveform in the linearized radiation theory 



is given by^ [43|, |4J] (see [ij, |45| for its derivation in the EFT framework) 



K, (*>x) = --A,,,,,^^-^^ + -^^n^ + -^^n„n^ + ...J (1) 

"' (^ 3 ^t^"" + 4 ^t^''"'''" + 15 dt^ ^aUmnn 

1 rl^ jl)bmnp 

where the time derivatives of the multipole moments are evaluated at retarded time tret = t — |x|, 
rii is the unit vector in the direction of observation and Aij^^i = PikPji — 2^'i-j-^ki is the transverse 
traceless projector with Pij = 6ij — nirij and («, j) implies symmetrization. 

As we discuss momentarily, this expression is sufficient for computing all of the spin effects in 



^ The overall minus with respect to the equivalent expression in [43l. |44| is due to the (+,—,—,—) metric signature 
we use. 



the waveform to 2PN and a sub-set of the 2.5PN contributions. At 2.5PN however, nonhnearities 



in the GW propagation first contribute to the spin components of t 
(scattering off the geometry) , previously computed in the hterature 
to this effect towards the end of the paper, and in more detail in 



le waveform via the tail effect 



e 37 
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491] . We will come back 



2.2. Multipole moments in terms of the stress energy pseudo-tensor 



The source multipole moments in the EFT formalism were given to all orders in terms of 
moments of the stress energy pseudo-tensor T^'^ in [14], where they were shown to coincide with the 
long wavelength expansion of the previously derived expressions for the source multipole moments 
in the more traditional FN formalism 46(]. For our purposes here the relevant results are 
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(2) 
(3) 
(4) 
(5) 
(6) 
(7) 



where the STF prescription only acts on the free indices of the multipole moments.^ Each moment 
of the stress energy pseudo-tensor is then computed from Feynman diagrams yielding T'^'^(t, k), 
defined as 



r^'^(t,k)= / (i'^xr'"^(t,x)e 



-ikx 



and using the Taylor expansion 13|, [35 1 



(8) 



T^^%t,]^) = Y, 



"'^" ' / d3xr^^(t,x)x*l • • -X^" ) k,, • • -ki^. 



n=0 



n\ 



(9) 



^ Note that the form of T^ differs from the expression in jl^], however they are physicaUy equivalent, which foUows 
from the conservation of yf". 
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b) 



FIG. 1: 0(Sa) worldline coupling to: a) radiation and b) potential modes. A blob represents couplings 
linear in spin. 



2.3. Overview of spin efTects in the EFT formalism 



In this paper we use the EFT formahsm for spinning compact objects developed in ll|, [l^, l22l . 
35]. The spin tensor, Sf^, is introduced in a local frame defined by the tetrad e^, which obeys 
^li^tVab = 9fiu and which is expanded in the weak field regime as e^ = J° + . . .. The spin vector is 
defined as usual: S"]^ = e"^^ S^. The coupling between gravity and spin reads [11] 



S., 



dXujfuf'Sl'', 



(10) 



where io?f are the Ricci coefficients and u^ = -^. In what follows we remove the subscript L for 
convenience. 

n 

For the radiation sector we follow our previous work [35], which we encourage the reader to 
consult for further details. The calculation of the spin contributions to the source multipoles is 
divided into two parts: worldline radiation and nonlinear contributions. The Feynman diagram 
for the radiation stemming solely from one worldline is depicted in Fig. [TJi. (Fig. ^ represents 
the coupling of a potential mode to a worldline 10, lUl] which only plays a role as part of a sub- 
diagram.) Expanding Eq. 
needed for our purposes is 



101), the interaction Lagrangian between spin and gravity to the order 

id, 
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22] 



^'^ = 2^^'^'^'"^'^'^" + 4^^"''^^ {l^-''^ + ^^ 



^^\,a-hf,a,x]v^' + 0{h•'), 



(11) 



Feynman rules 



11 



19|, 



d. 



where we choose the worldline parameterization X = t, such that u^ = v'^ = -^. This yields the 



351 ] needed for the spin couplings. 



The relevant expressions for the stress energy pseudo-tensor T'^'^(t, k) are [35 1 

T£{t,k) = Y,S''lik^e-''^-^^ (12) 

A 

^ft(*' ^^ = \Y. {^a i^' + S'i^A ^k^' + Sf zk • v^ - 5°^) e-"^-- (13) 

A 

W^ ^^ = \T.{ [Sa^a + S'}^\) ^k' + ^k . v^ {S^^^\ + S°V\) 
/I 

- S'^^\ - S'i^\ - S'^^\ - S'is.\] e-i^--, (14) 

where S^^'^ is the spin tensor that obeys the spin supplementarity condition (SSC) S^^ = kS'^^v^ + . . . 



11 



19l422l|. with K = 1 in the covariant SSC we adopt in this paper. Notice we suppressed the 
impUcit time dependence of the worldhne variables. Using Eq. ([9]), and expanding in powers of k 
to a given order, we read off the moments of the stress energy pseudo-tensor which enter in the 
multipoles in Eqs. (I2]17|). 

Nonlinear effects in the matching for the moments on the other hand are due to radiation off 
the gravitational part of the stress energy pseudo-tensor which accounts for the potential energy in 
the binary systems. The relevant diagrams thus involve couplings between radiation and potential 
modes as shown in Figs. [2}x&c. 
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TABLE I: Scalings for leading terms in moments of T^"^ for various orders in spin [35|. These scalings are 
valid for all moments £ > 2. 
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TABLE IL Order at which spin contributions to the multipole moments enter the GW amplitude. 

2.4. Power counting 

In order to know at which order spin effects start to contribute to the multipole moments and 
in turn to the waveform, it is instructive to derive the scaling of the leading contributions to the 
moments 



Kr 



d3xr^^(t,x)x^i...x^ 



(15) 



which we summarize in Table [J 35]. 

The scaling also clarifies which moments of T'^'^ in the expressions for the multipoles will con- 
tribute at a given order. From Table [I] and the expressions for the multipole moments^ we see that 
in order to compute the spin contributions to the waveform to 2.5PN we need, in addition to the 



results computed in 



only terms linear in spin. It is easy to show that 0{SaSb) terms do not 



enter until 3PN, and 0(8^) contributions are only needed to LO in the mass quadrupole, the mass 



octupole and the current quadrupole moments 



11 



35|. 



In [35| we computed up to NLO mass and current quadrupole and LO mass and current octupole 
moments linear in spin. Given the scaling shown in Table |T] we notice that current moments are 
enhanced by a power of v with respect to mass moments for terms linear in spin. For example, 
the LO and NLO /*■' enter at 1.5PN and 2.5PN in the expansion in Eq. ([1]), whereas J'-' appears 
already at IPN and 2PN respectively. For completeness we show in Table |II] at which orders the 



See [IJ, ISg] for the multipole moments beyond the ones in Eqs. (|2I7|I needed here. 



spin contributions to the multipole moments enter in the waveform. 

The remaining multipole moments we compute in this paper are all linear in spin, and consist of 
the LO J^^'''-, J«i'^'™ and P-'*^' as well as the NLO J*-'*^. Together with the other multipole moments 
in Table HI] already computed in [35.] and the tail contribution involving the leading J*-' linear in 
spin, described in Sec. 14.11 these complete the required ingredients for the GW amplitude to 2.5PN 
order. 

3. SPIN CONTRIBUTIONS TO THE SOURCE MULTIPOLE MOMENTS 

3.1. Current octupole J'^*^ 



Given its complexity we begin with the computation of the current octupole. As explained in 
35( we split the Feynman diagrams into worldline and nonlinear contributions. 



3.1.1. Worldline radiation 



We start computing the correction to Jq to NLO from the worldline couplings. Using Eqs. 
) and (fT3]) we find 



•^& = E 2 [SaX^a4] stf + f E [2viS^x^;,x^ - 2(Sa • ^a)^\^\A + ^^>\A{^A ■ xa) 



A A 

4S'^v^x^(xA • va) + (xA • SA)a'^x^x^ - {sla ■ xa)S^x^x^ 



STF 



(16) 



We now move to the expression for Jj-' in Eq. ([3]). The result reads 






(Sa • ^a)^\AA - 2xix>^S*;i + (xA • VA)Siix>^ - 3(vA • Sa)x^x^^x^ 



STF 

(17) 



Combining both, Eqs. ([T6]) and (fT7|) . we obtain 
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+2(xA • va)S^^v^^x^J ^ li " § J ^'^"^'aA^^A ■ XA 
+-(aA • Sa)x^x^;^x^ + -x^ [^\^\^\ + ^A^\^\ 



>^ ) ( v^x^;^x^(Sa • va) - v^S^x^;^xl - 2v^v^^x^(Sa • xa) 

K V 
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+ ^ )S^x^^x^(aA-XA) 



(18) 



STF 



for the worldline contribution to the current octupole to NLO. 



3.1.2. Nonlinear gravitational contributions 




m2 




FIG. 2: First nonlinear diagrams that contribute to J'^^^ at 0(8). As explained in the text, the diagram in 
Fig [5^1 represents a sub- leading effect. 



The relevant diagrams for the nonlinear gravitational contributions are depicted in Figs. [2jz6c. 
Since the expression in Eq. ([2]) entails the T'^* component of the stress energy pseudo-tensor, it is 
easy to show there is no contribution from Fig. [2]i at NLO. The same occurred in the computation 

r I I ' 

of the NLO J*-^ in 35(. Following the procedure described in ^] we now compute the remaining 
diagrams. The contribution from Fig. [2)o is relatively simple to calculate and reads 



GruB 
i^~ Z^ 27-3 

A,B 



^i = E ^ f(SA • XA)r'x^^x^ - (r • x^ + 8r2)S:,x^;,x^ 



STF 



(19) 



The most e 
is given by [35|] 



aborate calculation is that of Fig. [2]3, whose contribution to the stress energy tensor 



d^i 



l*^) - E S«X / ^i'P')i (k'k- + 2p'k. - *.'(k + p)^) ^J^e-'i^ 



hk)-xsgipx_4 



(20) 



In order to use Eq. (^ we need to compute J T^^x'^x^x^ , which fo 



by expanding to third order in k, as shown in Eq. ([9]) [13l . 
dropping irrelevant trace terms, we find 



lows from the previous expression 



35[. Performing the integrals, and 



■' A,B 



K^S^' U^'X + ^(x^^xl + xi^xl) + x^^xl 



(21) 



i^rr^r' Tx^^x^ + i(x>| + x^^x^) + x^^x^ ) + {m ^ j; m ^ k} - ^^^r^x^x^x^^ 
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Plugging this expression into Eq. ^ we obtain 
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J ^fc 



(SA-XB)rV^x^ + (r-XB)S^x],x^ + (rxSAr(xAXXB)^(x^ + x^) , (22) 



•?■ i'y'^. -U y1 
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which can also be re-written as 



jijk _ sr^ 
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i^j ^k 



STF 



(23) 



3.2. Mass and current 16-pole(s) 7*^'=', J'^'^' & current 32-pole J*^ 



fc^?^ 



The remaining multipole moments are needed at LO linear in spin and are given by worldline 
radiation alone. The computation follows the same steps as before, and the results are 



jijkl 



J^ 4 U + - (VA X Sa)*x^^x^x^ - — (XA X Sa) V^^x^xi 
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^dt ~ 2 ^ [^A^A^A^A 
A 

jijklm _ oST- \ci „i „A: I 
Odt ~ 2^ [^A^A^A^A^ 



STF 



(24) 
(25) 
(26) 



In the covariant SSC k = 1, these coincide with the ones given in Eq. (3.6) of [34i | for the appropriate 
order i in the multipole expansion. 
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TOWARDS SPIN EFFECTS IN THE WAVEFORMS TO 2.5PN ORDER 



Collecting all the ingredients, within the covariant SSC (k = 1), we finally obtain 

r^^ = Y, 2 [S^.^A^I] g^p + E \-l {^A^A^AiSA ■ v^) - viS^x^j^x^ - 2vi,V;,xl(S^ • x^) (27) 
A A ^ 

+ 2(x^ • Va)S\^'a^a) + Q^A^A^Ai^A ■ xa) - -S^x>^(aA • xa) 
+ g (aA • Sa)x^x^;^xJ^ + -xi (v^v^^Sl + x>^^S5i 
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-(r • XA - 8r2)S^x^^x^ + - (xiix^;ixl - x^x^^x|) Sa • r - -x^x^'^Sa • xa 



j.ifc^ = I E [4(VA X SaT^a^Wa - 3(XA X SA)V>ix!4 
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These source multipole moments, together with those reported in [35i] and including the necessary 
spin- independent ones, allow us to obtain the spin dependence in the gravitational waveform to 
2.5PN order from Eq. ([T]) using the equations of motion that derive from 



IC 
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Still 



missing is the contribution from the tail effect in the GW amplitude which we discuss next. 



4.1. Tail effect in the waveform 
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Let us now consider the corrections to the expression in Eq. ([T]) arising from nonlinearities 
in the GW propagation. In particular, we are interested in the tail effect which contributes to 
the waveform linear in spin at 2.5PN. The tail effect for GWs in the PN regime was first studied 



m 



47|, 



48( 1 where the tail corrections to the waveform were derived. In 



131 ] the tail effect was 



studied in the EFT formalism, however the results in [13] did not include the tail contribution 
to the waveform. We outline here how to extend the analysis in [13( to be applicable to the GW 
amplitude, and leave the more formal aspects for a future publication [421 ]. 



For our purposes we only need to consider the tail produced by the interaction between the 
current quadrupole J*-^ and the mass monopole M, shown in Fig. [3l This is because the tail 
contributes an extra factor of v^ compared to the leading linear contribution from the current 
quadrupole J*-' , cf. Table [III A similar diagram with J*-' — > P^ gives the leading order tail for 
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spinless objects at 1.5PN. 




FIG. 3: Tail diagram that describes the interaction between the source current quadrupole J*-' and the 
nionopole M . The dashed hne here represents the background mode which builds up the geometry sourced 
by the binary. 



Our starting point is the expression for the GW amplitude for a given source T*-^ 

hj^{t,^) = -^k,,M r^ ^T>'^{co,con)e'^'-^ + O (l/lxp) . (31) 

To obtain the waveform we need to compute the corresponding r*-'(a;,k = ojn), which is defined 
by means of the GW emission amplitude 

- ^e*^iun,h)T^^{u;,u;n) ^zAhiu^,un) = ^^, (32) 



Op 



} 5 



with ejj(a;n, ±2) the polarization tensor. For the case at hand, we need to compute the diagram 
in Fig. [3] in order to extract the relevant T*-' . Moments of the full stress energy pseudo-tensor 
computed this way, including tail effects, are often called radiative multipoles in the literature 44l |. 



Note that, formally speaking, the expression in Eq. ()3ip is only valid for sources with compact 
support. This is not the case for the stress energy pseudo-tensor induced by the J^^M interaction, 
which scales as 

r^'~l/|y|2 + 0(l/|y|3) (33) 

at a distance |y| away from the binary's worldline. However, the GW scattering off the long range 
Newtonian potential induced by the binary at distances far away from the source only leads to a 
shift in the phase which can be absorbed into an unobservable time redefinition, as will be seen 
below. For LIGO/Virgo, which are only sensitive to the last few minutes of the binary inspiral, 
the waveform must be first observed as it enters the frequency band to start tracking its evolution, 



13 



and it is this change in time which ultimately corresponds to a measurable quantity. This justifies 
the use of the expansion in Eq. (i3T]) . See [42] for a detailed discussion. 



The tail efi'ect in the graviton emission amplitude for the P^ M interaction was studied in detail 
in [13|] , represented by the diagram in Fig. [3] with J*-' -^ P^ . As shown in [13] the corresponding 
integral is IR divergent and can be regularized using dimensional regularization. The result reads 



i Ajij ^j {en , con) = iAji] x (iGMca) 



juj + ieY' 



alE 



(d-4)/2 



d-A 
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(34) 



where Ajij = j^io'^e*jP^{io), and je = 0.577... is the Euler constant. Note there is a slight 



difference from the result in [l3|], i.e. the factor of (a; + ie)^ versus cj^ + ie in Eq. (18) of 13|. 
because of the necessary retarded boundary condition for the propagator, which takes the form 
i/ [{uj + ie)^ — k^) in Fourier space. 

From here we can read off the contribution from the tail to the radiative quadrupole moment 
by expanding around d = A (we absorb the factor of vr into /i from now on), 



^Idi"^) = ^o''(^) { 1 + GMcj ( sign(a;)^ + i 



2 uj^ 
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h log ^- + 7s - 




eiR IJ'^ 
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(35) 



Notice the presence of an infinite IR divergent piece proportional to 1/eiR with e = d — 4 and the 
presence of the arbitrary scale /i. The IR divergence arises due to the effect of the Newtonian 1/r 
potential sourced by the total (ADM) mass of the binary on the propagation of the GW. Similar IR 
divergences also appear in electrodynamics (QED). As is well known, the 1/r Coulomb potential 
produces an infinite phase shift in scattering amplitudes [5D:]. However, this is unphysical in the 



sense that we do not observe the absolute va 
out of any measurable quantity. See also [5C 



ue of the phase shift in QED, and the infinities cancel 



511 ] for similar IR effects in gravitational scattering. 



As any IR divergence, the one which occurred here in the computation of the tail cannot contribute 
to physical quantities. In fact, one can show the IR divergences exponentiate into a pure phase 
in the emission amplitude, i.e. A = ^iR-finite e^**^^^'^' "^'^ , and therefore drop out of observables 
extracted from ].4]^, such as the total power [131 ]. 

In the waveform on the other hand, the l/ejR pole remains as an infinite overall phase. However, 
it is easy to show we can absorb it into a time redefinition, namely 



tha 



i-2GM/eiK, 



(36) 
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where tbare is the (infinite) bare time, whereas t represents the time variable used by the exper- 
imentahsts. This is in accord with our previous discussion, in the sense that one must choose a 
time origin to start tracking the signal, which is of course arbitrary. As we will show in detail in 
42l ]. the IR divergent contribution can always be absorbed into a shift in time. 

This renormalization procedure is arbitrary, and that is the reason for the fi scale in the above 
expressions. Furthermore, one can easily show that a change fiQ — )• A^o can be cancelled out by a 
concurrent transformation t(A/io) -^ ^(/^o) + 2GMlog A. Hence the GW amplitude is independent 
of fl. 

From the above expressions the contribution from the P^ M coupling to the GW amplitude 
becomes 



{k 



ij )pJ+IiJM ^*'-'^'' 



2G, 

■-p-rA,- 
X 



ij,kl 



^^ 27r 



2gia;t,et(/.)gi0taii(<^,M) (1 + GM\uj\-k)\ rj{uj), (37) 



-ijj e 



with 



^2 -^Y 

7taii(w, /u) = GMu) ( log — + 7£; - — ) . 



(38) 



Up to an unphysical rescaling of ;U, this agrees with the result previously obtained in |48]. A sim- 
ilar computation gives us the contribution from the current quadrupole and monopole interaction 
which leads to the leading order spin-dependent effect from Fig. [3l and reads 



( 



^ij IjiJ+JiJM^'^'^' 



2G . 

X 



^^ duj 
^ 2^ 



_^2gia;t,et(M)gi^taii(<^,M) (1 + GMlwIvr)] fj{u), (39) 



with 



^2 J 

V'taiKw, ^J) = GMu ( log — + 7B - - 



(40) 



This agrees also with the result in [37|, |4£ 

Note that while the constants associated with the logarithms in (j38p and (j40p are unphysical, 
the difference is not, since we are only free to make one choice for /x. The 7£;'s are universal (i.e. 
independent of the multipole) and thus unphysical, but the rational number is non-universal as 



can be seen from (1381) and (HOl). see also 



44| 



The independence of Eq. (|39p from any particular choice for fi guarantees physical observables 
do not depend on any specific set of coordinates, as one would expect. However one still has to 
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m, 



make a choice, for example fi = uig (the seismic cutoff frequency for earth-based detectors 
which will be directly linked to a reference frequency scale at the time to the GW signal is first 
observed as it enters the detector's band. Once a choice for fj, is made, and the initial value of the 
GW phase is measured, the expression in Eq. ([391) tells us how to correlate the predicted evolution 
of the GW with the observed signal. We will formalize these arguments in greater detail in a future 
publication |43]. 

5. CONCLUSIONS 



In this paper we computed the multipole moments, expressed in Eqs. (1271- 1301) in the covariant 

y, y, 



m, 



mm, 



35l | lead to the spin 



SSC, which together with the results reported in 

contributions in the GW amplitude to 2.5PN order. We employed the EFT formalism developed 

in 10l4l4l] , with which we previously computed all the spin-dependent ingredients that enter in the 



GW phase to 3PN order in [35]. We also reproduced the leading order spin contribution to the 



tail effect, already obtained in the literature J37l. l47l| . In a separate publication [53(, we will report 
the expressions for the GW polarizations with spin effects to 2.5PN order in a ready-to-use form 
amenable to GW data analysis. 
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